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ANALYSIS OF A HYBRID PHASE-LOCK LOOP 


By Thomas L. Stewart 
Electronics Research Center 


SUMMARY 

A model for a hybrid phase-lock loop has been derived by 
using Fokker-Planck techniques. The hybrid loop operates on both 
the data and carrier components of the received signal to provide 
the phase estimate required for coherent detection. The loop is 
optimized with respect to the system design parameters, and the 
performance of a bipolar phase-shift-keyed communication system 
is evaluated. The results are compared with systems which employ 
Costas-type loops or conventional phase-lock loops. The hybrid 
loop gives superior performance for all signal-to-noise ratios, 
with a maximum improvement of 1.5 dB for low data-rate systems, 
and, additionally, is satisfactory over a wider range of signal- 
to-noise ratios than either the Costas or conventional phase- 
lock loops. 


INTRODUCTION 

Binary phase-shift-keyed (PSK) communication systems transmit 
information in the form s (t) = /2 1 A cos (wt+a (t ) 0 o +0 ) , a(t) = ±1. 

To perform coherent detection, it is required that the receiver 
determine the carrier phase 0 prior to detection. If the signal 
s (t) contains a residual component at the carrier frequency, a 
phase-lock loop can be used to track the carrier phase. On the 
other hand, the modulated portion of the signal also contains 
information about the carrier phase ; systems for generating a 
phase reference directly from the modulation component have been 
described by Costas (ref. 1) . 

This report is concerned with a phase-estimation scheme 
that combines the phase estimate from the carrier component with 
the phase estimate from the modulation component. It is mechanized 
by employing a combination of the conventional phase-lock loop 
and the Costas-type phase-lock loop. This combination is called 
a hybrid phase-lock loop (hereinafter referred to as hybrid loop) . 
The remaining sections of this report are devoted to a description 
and theoretical analysis of the hybrid loop. Fokker-Planck tech- 
niques are applied to measure its performance and the results are 
used to evaluate the error probability of a binary, bipolar 
communication system. Detection probabilities for the hybrid loop 
are compared with those for equivalent systems which use a phase- 
lock loop or a Costas-type loop. 
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THE HYBRID-LOOP 

The received signal is assumed to be of the form 
y (t) = s (t) + n (t) (1) 

with 

S(t) = /2A sin (cot + x(t) cos ^ m + 0) (2) 

and 

n(t) = n (t) cos cot + n (t) sin cot (3) 

c s 

where A 2 is the received signal power; x(t), a bipolar modulated 
square-wave subcarrier of unit magnitude; m, the modulation index 
which adjusts the power allocation between the modulated and 
carrier portions of the total signal; 0, the phase shift intro- 
duced during transmission and assumed unknown at the receiver; 
and n(t), the additive noise component. 

Expanding s (t) gives 

s(t) = /2Am sin (cot+0) + /2A (1-m 2 ) x(t) cos (cot+0). (4) 

From Eq. (4) we note that the data component with power A^l-m^) 
is 

s^(t) = /2A (1-m 2 ) 1//2 x(t) cos (wt+0) (5) 

and the residual carrier component, which will be termed the 
reference signal, is 

s r (t) = /3Am sin (ojt+0). (6) 

The physical significance of the parameter m may be seen from 
Eq. (6) where m may be defined as the square root of the ratio 
of the power in the reference signal to the total transmitted 
power, i.e., m = (P r /P) 1 /2. it f s worth a digression at this 
point to briefly review several schemes for generating a phase 
reference at the receiver that will lead us to the technique 
under consideration here. 

For the case in which m > 0/ Viterbi (ref. 2) has performed 
an analysis by Fokker-Planck techniques of a scheme whereby the 
carrier or reference component can be tracked by a phase-lock 
loop. An experimental evaluation of this method has been made 
by Lindsey and Charles (ref. 3). Stiffler (ref. 4) and Lindsey 
(ref. 5) , moreover, have considered the question of optimizing 
the modulation index m to minimize total system-error probability. 
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For m = 0, in the absence of a reference component, Costas and 
squaring-loop mechanizations (refs. 1 and 6, respectively) have 
been proposed. Didday and Lindsey (ref. 6) have demonstrated the 
mathematical equivalence of these two systems and have applied 
Fokker-Planck methods to analyze their behavior. Bussgang and 
Leiter (ref. 7) have derived the performance of a communication 
system, which uses a reference signal and an optimum detection 
scheme. This is equivalent to employing a maximum likelihood 
estimate of the phase in a correlator-type detector. Finally, 

Van Trees (ref. 8) has considered a system, based on linear track- 
ing theory, that is similar to the one proposed here but is re- 
stricted to high signal-to-noise- ratios to permit linearization 
of the loop equations . 

The hybrid phase-lock loop is shown in Figure 1. As stated 
in the Introduction, it is a combination of the Costas and con- 
ventional phase-lock loops. By setting the gain a 2 to zero, the 
hybrid loop reduces to a conventional phase-lock loop which de- 
rives the phase estimate from the carrier component of y(t). 
Setting ai to zero gives a Costas-loop mechanization which ex- 
tracts the phase from the data component of y(t). In essence. 


MULTIPLIER 

Ml 



M2 


multiplier 

M3 


a,, a 2 = LOOP -GAIN FACTORS 

H,(s) f H 2 (s), H 3 (s) = LOOP 'FILTER TRANSFER FUNCTIONS 
Figure 1.- Block diagram of hybrid loop 
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the hybrid loop combines the error signals from both loop branches 
to drive the voltage-controlled oscillator (VCO) . The loop fil- 
ters with transfer functions H]_(s) , H 2 (s) , and H 3 (s) remove cer- 
tain double frequency and modulation components resulting from 
the multiplication operations, and are otherwise assumed to be 
flat over the frequency range of the loop. 

In addition to the usual system design parameters, the hybrid 
loop must be optimized with respect to a^ and a£ , the two gain 
factors which provide the relative weighting of the hybrid-loop 
error signals. Also, the modulation index m must be determined 
to permit proper allocation of the transmitter power between the 
data and carrier components . 


ANALYSIS OF THE HYBRID LOOP 


In the configuration depicted in Figure 1, the outputs of 
multipliers Ml and M2 after filtering by Hi(s) and H 2 (s) are 
given by 


e o (t) 


-C'-' - n s (t) , - 

Am sin (0-0) + — cos 0 - sin 0, 


n c (t) 




e 1 (t) 


2 1/2 - n c (t) 

All-mT 7 x ( t) cos (0-0) + — cos 0 




n (t ) 


s 


VT 


sin 0 , 


and 


e 2 (t) 


A(l-m 2 ) 1/2 x ( t ) 


A 

sin (0-0) 


n c ( t) 


S\ 

sin 0 


cos 0, (7) 


where it is assumed that Hj(s) and Hjfs) remove all double-fre- 
quency components due to the multiplication. It is further 
assumed that the bandwidth of Hi (s) is sufficiently narrow to re- 
move the modulated subcarrier x(t) and that H 2 (s) removes the 
low-frequency term from the carrier component. According to these 
assumptions it would be required that the data subcarrier x(t) 
have insignificant spectral components near the origin. 


n s (t ) 

V3 
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The output of multiplier M3 is given by 


e 3 (t) 


A^l-m 2 ) 1 / 2 

2 


n„ (t) 


sm 


c vw 2 1/2 

2<j> + A(l-m ) ' x(t) sin (<j>-0) 

V 2 " 


n s 2 1/2 ~ 

- A(l-in ) x/ x(t) cos (<f>-0) 


n s (t)n c (t) 


’n 2 (t) n 2 (tf| 

4 4 j 


sin 20 


cos 2 0 


( 8 ) 


where <j> = 0 - 0 . The VCO error voltage is 


e(t) = a 1 e Q (t) + a 2 e 3 (t) 


(9) 


where the factors a]_ and a 2 combine the estimates from the two 
sections of the loop. Let us denote the VCO constant by K, which 
also includes any additional loop gains; then the loop differen- 
tial equation may be written 


d(j> 

dt 


-K e ( t ) 


( 10 ) 


in which the detuning of the VCO from the carrier frequency to is 
assumed negligible. 

Equation (10) has the form cj) = F(<J>,t), and the corresponding 
Fokker-Planck equation is, in the stationary case. 


1 

2 



k 2 (c!>) p (< p ) 


it P ( « 


0 . 


( 11 ) 


Here 


K 1 (4>) = E F ( <j) , t ) , 


(12) 


5 



(13) 


ym OO 

K 2 ( < t>) = J E F (<j) , t) F (<t>,t+T) 


K 2 dx. 


E denotes ensemble averaging, and the assumption is made that the 
correlation time x c of the noise process is much smaller than the 
loop time constants.* Physically, this is equivalent to the re- 
quirement that the loop bandwidth is much smaller than the band- 
width of the noise in the loop. Details of this concept are 
given in reference 9. 

If we now assume that cf) (t) and x(t) are slowly varying so 
that x(t) = x ( t+T ) , then 


K x (4>) = E F (<p , t) = -F Ke (t) 

a^Am sin cj) + a 2 2 


= -K 


A 2 (1-m 2 ) 


sin 2(p 


] 


(14) 


and 


K, 


= f 4 f l R nn< T > + if [' 

•— OO V- 


-A 


R nn (T) + "I | R nn (T)+2A2(1 - m2)R nn (T) 


’dx . 


(15) 


Substituting Eqs . (14) and (15) into Eq. (11) and solving for 

p (4>) , we have 


p (<t>) = C exp 


2K 

K, 


2 2 

a-^Am cos cj) + a 2 — — ^ m - cos 2<p 


(16) 


where C is determined by the condition 


The correlation time of the random process n(t) is defined as 


c 

t c 


R nn (T) 

R nn (0) 


dx 


where R (x) is the correlation function of n(t). 
nn 
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(17) 



d<j> 


1 


Next, let the envelope of the input noise process have a 
flat power-spectrum density with magnitude N 0 W/Hz and cutoff 
frequency to. Hz. The noise correlation function is 


R n n <T > = R n n < T > = 2N o“i 
c c s s 


sin 2 tho^t 
2ttuk t 


(18) 


and substitution of Eq . (18) into Eq. (15) gives 


K 2 N 

K 2 = 2 


T- 2 [ a J + a 2 N o"i + a 2 A2(1 - ”> 2 )] 


(19) 


From the linear phase-lock loop theory, the equivalent noise band- 
width of the loop is defined as 


2 or 


= / 


H (s) ds. 


( 20 ) 


-D‘ 


The expression for H(s) is found by linearizing Eq. (10) giving 


H ( s ) = 




KA| ( 1-m ) 


] 


s + KA.|^a^m+a 2 A ( 1-m 2 )J 


( 21 ) 


which yields 




ka r 
4 [ a l 


m+a 2 A(l-m )J Hz. 


( 22 ) 


With the substitution from Eqs . (20) and (22) into Eq. (16), we 

have 


{• 


p (<f>) = C exp cos <f> + cos 2c}> 




(23) 
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where 


D = a6 [l+pyjm'' 


(24) 


and 


E = 1 + p‘ 


r. — h>] 

L a 6 3 ( 1 -ni )J 


(l-m z ] , 


(25) 


in which expressions 


A* 

a NR 
o 


6 = 


0 ), 


a 2 A 


Y = 


l-m‘ 

m 


Ka,Am u- 

“ L = -?— u+pyI e = -. 


(26) 


With set equal to zero, or p = 0 , Eq . (23) reduces to 


p(4>) 


C exp a6m cos 4> , 


(27) 


which is the phase-error probability density function for a con- 
ventional phase-lock loop as described in references 2 and 3. If 
a^ is set equal to zero, or p = 00 , then Eq. (2 3) reduces to 


p ( cj> ) = C exp D cos 2c(), 


(28) 


with 


D 



aS 



(29) 


which corresponds to the probability density function for the 
Costas loop (ref. 6) . 
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PERFORMANCE OF THE HYBRID LOOP 

In order to gain some insight into the degree of improvement 
offered by the hybrid loop, we will first compare the variance of 
the hybrid-loop phase error with the variance of the phase error 
for the Costas and conventional phase-lock loops. This will be 
done under the assumption that the probability density function 
of the phase error can be approximated by the Gaussian density 
function. Then, we will obtain exact results for the performance 
of correlation detectors. 

Approximating cos <j) by(l - cf> 2 )/2 in Eqs . (23), (27), and (28), 

we find the variance of the three loops to be as follows : 

1 . Hybrid Loop 

1 + p 2 (1-m 2 ) 1 + ~- 

2 L a<5 8 (1-m ) 

CT H „ 2 

aS [m+p ( 1-m ) ] 


(30) 


2 . Conventional Phase-lock Loop 


2 _ 1 _ 

^P 2 

ouSni 


(31) 


3 . Costas Loop 


2 

CT C = 


1 + 


a6B (1-m ) 




aS (1-m 2 ) 


(32) 


The optimum value of p for combining the loop outputs can be 
found from Eq. (30) to be: 


a6 3 (1-m ) 


opt m |l+a6B (l-m 2 )J 


, 2 
P . 1-m 

■? m 


(33) 
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MEAN -SQUARE PHASE ERROR 


Thus, for fixed m, the different loop outputs are combined accord- 
ing to the ratio of the phase-error variances from each section of 
the loop. If Eq. (33) is substituted for p in Eq. (30) , we have 
for the hybrid loop 


2 _ _1 [ l+a5 g (1-m 2 ) ] 

H a6 f 2 2 2 2 ' 

Ia6(3 (1-iu ) + irr[l+a<5g (1-nT) ] 


From Eq. (34) , we see that the variance of the hybrid loop is 
always less than that of the Costas or phase-lock loop. Also, 
the maximum improvement is a factor of 2, which occurs when 

CTp = O c* 

The variance of the phase error is plotted in Figures 2, 3, 
and 4 for various combinations of the loop parameters. Examina- 
tion of the curves indicates that the performance of the phase- 
lock loop is superior at low signal-to-noise ratios, while that 


4 + 4 


(34) 



LOOP SIGNAL-TO-NOISE RATIO, a8,dB 


Figure 2.- Mean-square phase 

error a 2 versus loop signal 

to-noise ratio aS = A 2 /n 0 wl 

for hybrid, Costas, and con 

ventional phase-lock loops 

with B = m I to . = 0.1 
L 1 i 
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MEAN -SQUARE PHASE - ERROR , tr 2 MEAN -SQUARE PHASE ERROR, 



LOOP SIGNAL -TO -NOISE RATIO, 08, dB 


Figure 3.- Mean-square phase 
error o 2 versus loop signal- 
to-noise ratio a& = A 2 /N 0 col 
for hybrid, Costas, and con- 
ventional phase-lock loops 
with 8 = to I to . =0.05 

-Li 1 



LOOP SIGNAL -TO -NOISE RATIO, a 8, dB 



Figure 4.- Mean-square phase 
error a 2 versus loop signal- 
to-noise ratio a6 = a2/n o 0)l 
for hybrid, Costas, and con- 
ventional phase-lock loops 

with 3 = oj t I a) . = 0.01 
L 1 i 
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of the Costas loop is better at high signal-to-noise ratios. The 
performance of the hybrid loop is found to be superior to either 
the Costas or conventional phase-lock loop, with the maximum im- 
provement at the point where the plots of and intersect. 

Of practical importance is the fact that the performance of the 
hybrid loop is better over the entire range of signal-to-noise 
ratios in contrast to the Costas loop, which is poor at low sig- 
nal-to-noise ratios, and the phase-lock loop, which is poor at 
high signal-to-noise ratios. Consequently, the hybrid loop will 
be less sensitive to variations in the signal-to-noise ratio over 
its operating range. 

We now want to compute the performance of a correlation 
detector which uses a hybrid loop for the phase estimate. The 
correlation detector computes 

/* T 

A = / y(t) x(t) cos (cot+0 ) dt . (35) 

0 


If A > 0, the decision is made that the binary data is +1, and -1 
if A < 0. It is easy to show that the detector-error probability 
P(e), conditioned on (p , is 

l ~x 2 /2 

— e dx . (36) 

cos <j> 


The average error probability is found by averaging P(e|4>) over 
p(4>) given by Eq. (23). That is, 


P(e(cf>) = 


/ 


< 


— 2 2 “ 

2A Z (1-irT) 

N R 
o 


P (e) 


/ 7T 

PU| 


<f>) p (<f>) dcj> . 


(37) 


The system is optimized with respect to m and p, where p opt 
and ni 0 p t are those values of m and p which minimize P(e) . A 
computer program was written to determine p opt and m opt • The re_ 
sults are presented in Figures 5, 6, and 7. ^Figures B, 9, and 10 
show the corresponding error probabilities. 
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OPTIMUM VALUE OF m OPTIMUM VALUE OF 




Figure 5 . - Optimum parameter 
p and optimum modulation index 
m versus input signal-to-noise 
ratio a = A^/N 0 R for hybrid 
loop with 6 = a) L |to. = 0.1 



0.6 P 



INPUT SIGNAL-TO-NOISE RATIO, a=A*/N 0 R,dB 


Figure 6 . - Optimum parameter 
p and optimum modulation index 
m versus input signal-to-noise 
ratio a = a2/n q R for hybrid 

loop with B = or. I to. = 0.05 
L 1 i 
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R/tu L = 12 
• 20 
= 28 



0 12 3 4 

INPUT SIGN AL-T0- NOISE RATIO, a =A 2 /N 0 R , dB 


Figure 7.- Optimum parameter 
p and optimum modulation index 
m versus input signal-to-noise 
ratio a = A-2/N 0 R for hybrid 
loop with 3 = 1 60 = 0*01 



-3 -2 -10 123456789 

INPUT SIGNAL-TO-NOISE RATIO, o=A*/N 0 R,dB 


Figure 8.- Receiver error prob 
ability P(e) versus input 
signal-to-noise ratio 
a = a2/n q R for hybrid loop 
with 8 = co L | = 0.1 



ERROR PROBABILITY, P(c) E „ RR0R PB0BABIL,TY . 



INPUT SIGNAL-TO- NOISE RATIO, a =A f /N 0 R, dB 


Figure 9.- Receiver error prob- 
ability P(e) versus input 
signal-to-noise ratio 
a = A^/NqR for hybrid loop 

with g = to T I co . =0.05 
L 1 l 



-2-10 I 23456789 
INPUT SIGNAL-TO-NOISE RATIO, a =A*/N 0 R, dB 


Figure 10.- Receiver error prob 

ability P(e) versus input 

signal-to-noise ratio 

a = a2/n d R for hybrid loop 

with 6 = a>_ I ai . =0.01 
L 1 i 


In Figures 11, 12, and 13, the hybrid loop is compared with 
equivalent systems which use a Costas or conventional phase-lock 
loop. The error probability for the phase-lock loop system is 
given by Eq. (37), with p(<j>) given by Eq. (27). The phase-lock 
loop system is optimized with respect to m. For the Costas loop, 
the phase estimate exhibits a phase ambiguity of tt radians. Con- 
sequently, a realistic comparison of this system must include the 
additional degradation due to this difficulty. One scheme for 
resolving the problem employs differential encoding of the trans- 
mitted data (ref. 2). The equation giving the error probabil- 
ity for this system is 


tt/2 

P(e) = f 2 P (e | cj)) [1-P (e | <f>) ] p(4>)d<f> 

- tt/2 

with p (4>) as defined by Eq. (28), and P(e|<t>) by Eq. (36) 
m = 0 . 


(38) 
wi th 



Figure 11.- Comparison of error 
probabilities versus signal-to- 
noise ratio for hybrid, Costas, 
and conventional phase-lock 
loops with 8 = w T jo). = 0.05 and 
6 = R/w l =12 L 1 
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• HYBRID LOOP 0 = 0.01 
m = rT1 opt P S P opt 
- COSTAS LOOP0 =0.01 
m = I Differential Encoded 
-PHASE-LOCK LOOP 


Figure 12.- Comparison of error 
probabilities versus signal-to- 
noise ratio for hybrid, Costas, 
and conventional phase- lock 
loops with 3 = to loo. = 0.01 and 
6 = R/w l =12 L 1 


-2 -10 1 234 5678 

INPUT SIGNAL -TO - NOISE RATIO, a=AVN 0 R, <*B 


- HYBRID LOOP 0 = 0.01 
P~P opt m 5 m opi 

- COSTAS LOOP 0 = 0.01 
m=l Differential Encodin 

- PHASE-LOCK LOOP 
m = m opt 


Figure 13.- Comparison of error 
probabilities versus signal-to- 
noise ratio for hybrid, Costas, 
and conventional phase-lock 
loops with 3 = oo | w. = 0.01 and 
6 = R/ co =28 L 1 

±J 


-2-10 I 23456789 
INPUT SIGNAL-T0 - NOISE RATIO, a = A f /N 0 R, dB 



CONCLUSION 


The hybrid loop offers significant improvements when compared 
with equivalent systems utilizing a phase-lock loop or Costas 
loop. The first improvement is that for low data-rate systems, 
R/wl < 20, the hybrid-loop performance is less sensitive to varia- 
tions in the signal- to-noise ratio. The same is not true for the 
Costas loop, which degrades rapidly for signal-to-noise ratios 
less than 3 dB, or the conventional phase-lock loop, which is in- 
ferior for signal-to-noise ratios greater than 3 dB. The hybrid 
loop is, in fact, uniformly better than either of these methods 
and would give additional protection to systems which operate 
under wide variations of the signal-to-noise ratio. 

In addition, the hybrid loop uses the total received signal 
power to provide a phase estimate. Depending on the operating 
signal-to-noise ratio, the hybrid loop can offer a gain improve- 
ment in performance of about 1.5 dB when compared with systems 
which employ a phase-lock loop or Costas loop. 
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TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 
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